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I. INTRODUCTION 

As a natural extension of early work on relativistic fluid mechanics 1 , Lichnerowicz and 
Anile developed a theory of relativistic magnetohydrodynamics 2-4 paralleling the well- 
studied nonrelativistic version. The primary assumption of MHD, that the fluid in question 
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is charged but quasineutral, holds in relativistic contexts of interest, although the definition 
of quasineutrality must be restated in terms of a 4-current. As a result, relativistic MHD 
holds an important position in the field of relativistic computational modelling, with a 
variety of algorithms both suggested and implemented (e.g. Refs. 5-11). The present paper 
explores the theoretical side of the subject, which recently has received less attention than 
the computational side. In particular, our main contributions are i) the introduction of a 
new canonical 4-momentum, and a new divergenceless 4-vector to represent the magnetic 
field; ii) using the new variables to cast relativistic MHD into a covariant Poisson bracket 
formalism in terms of Eulerian field variables; iii) investigating many properties of our new 
formalism, including several alternative brackets, a reformulation in differential-geometric 
concepts, and the consequences of a new gauge freedom. 

Physicists know well the benefits of casting a theory into a Hamiltonian or action prin¬ 
ciple mold, as our present work accomplishes. In addition to being aesthetically appealing 
in its own right, this form has several practical advantages: (i) certain numerical algorithms 
are based on such a structure (e.g. the recent works of Refs. 12-14), while others can use said 
structure as a consistency check; (ii) finding the equations of motion in general coordinates, 
which Landau and Lifschitz called “unsolved” for fluid mechanics, becomes straightforward; 
(iii) the formulation assists both the discovery and classification of constants of motion; (iv) 
a Hamiltonian structure provides a handy framework for equilibrium and stability analysis; 
(v) both Hamiltonian and action principle pictures provide a way of quantizing physical 
systems, tying into the field of “quantum plasmas” currently receiving much attention. The 
present work also necessitates a handful of new concepts (a modified enthalpy density, a 
momentum differing significantly from the standard kinetic momentum, and another “mo¬ 
mentum” conjugate to the magnetic field), which may provide new insight into this physical 
system. 

It would be remiss for us not to mention previous attempts at providing an action principle 
for relativistic MHD. Maugin 15 did provide a Lagrangian action principle, but in terms of 
Clebsch potentials, rather than the physical quantities themselves. Alternatively, Kawazura 
et al. 16 have recently produced a useful Lagrangian variable action principle. In future 
work we will show how these action principles relate to our own. Meanwhile, the Poisson 
bracket structure of Morrison and Greene 17 was shown to be applicable to relativistic MHD 
in terms of 3-vectorial quantities in a specific reference frame 18 . This bracket is in effect a 
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(3+1) split of the present theory, which uses only tensorial quantities and does not require 
a choice of reference frame. The chief advantage of the present work over Maugin’s and 
Kawazura’s is that it takes place in Eulerian variables, rather than Lagrangian ones: both the 
aforementioned theories require a map to the physical Eulerian variables after the variational 
principle has been performed, adding an additional step that our formalism does not require. 

The paper is organized as follows. Section II provides a review of MHD, starting with the 
nonrelativistic theory in Sec. IIA before describing the relativistic theory in Sec. IIB, where 
the new variable /+ that describes the magnetic field is introduced. Section III then presents 
our new action principle using the new variable. Here we first describe, in Sec. Ill A, the 
functional that serves as our action and show how conjugate variables arise from functional 
differentiation; then, in Sec. Ill B, we describe the covariant Poisson bracket formulation 
that provides our constrained variations. Section IV is dedicated to alternative brackets: 
first, in Sec. IV A we present a bracket with nontrivial Jacobi identity; in Sec. IV B, a bracket 
using a tensorial magnetic potential; in Sec. IV C, the differential-geometric form of our main 
bracket, with several other quantities also presented in that form; finally, in Sec. IV D, we 
show how to couple our relativistic MHD theory to a fixed gravitational background. In 
Sec. V we discuss several features of our theory, including the nature of the divergence-free 
constraint and Casimirs. Finally, in Sec. VI we summarize our results. 


II. MHD REVIEW 


The equations of MHD, both nonrelativistic and relativistic, can be written in various 
ways in terms of different variables. In this section we gather together formulae and well- 
known identities needed for the remainder of the paper. The main new contribution of this 
section is the introduction of the variable hJ J of (12). 
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A. Nonrelativistic MHD — two descriptions 


First we give the equations of ideal nonrelativisitic ideal MHD, with the force law and 
Faraday’s law expressed in two alternative ways: 


<9v 

dt 


(v ■ V) v = [ (V x B) x B] 


p An p 

Vp 1 

- 1 - — 

p Anp 


V- (IB 2 / 2-B<8)B) 


HD 

_ = V X(VX B ) 

= -B V V + B Vv - V VB 

| + V(pv)=0 

ds ___ „ 

_ + V • Vs = 0. 

at 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


Here p is the fluid density, p its pressure, s its specific entropy, v the velocity field, and B 
the magnetic field. In (2) the symbol Z represents the identity tensor. The current j and 
electric held E have been eliminated from these equations, but they can be recovered from 
the ideal conductor Ohm’s Law, E + (v/c) x B = 0, and Ampere’s Law, j = (c/47r)V x B. 

Observe the alternative versions of (1) and (3) given in (2) and (4), respectively. These 
equations differ by terms involving V • B, and both Eqs. (3) and (4) preserve the initial 
condition V • B = 0, which can be seen by rewriting (4): 

— = -B V ■ v + B ■ Vv - v ■ VB = V x (v x B) - v V ■ B . (5) 


Upon taking the divergence, 

^ = -V-(W-B). (6) 

Consequently, if V • B is initially identically zero it remains so as well. Equation (5) shows 
that forms (3) and (4) are equivalent when the magnetic held is divergenccless, although 
the former reveals its Faraday law origin, while the latter shows an advected magnetic hux. 
Geometrically (4) is <9B/<9f+X v B = 0, where X V B is the Lie derivative of B, a vector density 
dual to a 2-form. Similarly, Eqs. (1) and (2) differ by a V • B term, with the former revealing 
its Lorentz force origin via a clearly identihed current, while the latter takes the form of a 
conservation law, which Godunov 19 showed to be superior for numerical computation. 
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We have distinguished these two forms because they possess different Hamiltonian struc¬ 
tures. In Ref. 17 a Poisson bracket was given for the form with (1) and (3), but this structure 
required building in the initial condition V • B = 0. However, an alternative and more nat¬ 
ural form was first given in Refs. 20 and 21, which is entirely free from V ■ B = 0, it being 
only one possible choice for an initial condition. Later in the paper we will demonstrate 
relativistic equivalents of both structures, and the two will also differ by the divergence of a 
4-vectorial quantity; to be equivalent, said divergence must vanish, which will motivate our 
use of the new magnetic quantity h A 

Should one wish to add displacement current back into MHD, as is done in the most 
prevalent version of relativistic MHD, the momentum equation would have to be altered as 
follows: 


<9v Vp 1 

— + v ■ V v =-- + — 

at p 47t p 


V x B + 


d_ 

dt 


x B x B 


(7) 


However, the new term, when compared to dv/dt, scales as 


B 2 


An pc 2 


va 


where Va is the Alfven velocity. In the nonrelativistic limit, waves involving disturbances of 
the matter must also travel much slower than the speed of light, allowing one to drop the 
displacement current. This also means that relativistic MHD is free to add said displacement 
current back in (albeit constrained by Ohm’s Law), while still reducing to conventional MHD 
in the nonrelativistic limit: one simply needs to keep in mind that said limit goes beyond 
just setting v/c —> 0. 


B. Relativistic MHD 

Turning now to the description of relativistic MHD, we use signature and units such that 
4-velocities have positive unit norms = g^u^y? = 1, where the Minkowski metric g^ u is 
given by dia(1, — 1, — 1, — 1). The 4-vector field u M will denote the plasma’s 4-velocity at each 
point in spacetime; at each such point, this quantity will define a reference frame with locally 
vanishing 3-velocity, helpful for some purposes. The fluid density is now p = mn(l + e), 
where n is the baryon number density, m is the fluid rest mass per baryon (including both 
proton and electron, for the typical case), and e is the internal energy per baryon, normalized 
to m. The specific entropy s is unchanged, though later on it will prove more convenient 
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to use the entropy density a = ns. We will suppose that the energy can be written e(n, a), 
hence p(n, a), in which case the pressure is given by 

d P , dp 

p = n a^ + a B^- p ’ (8) 

which is just the first, law of thermodynamics, Tds = d(p/n ) +pd(l/n), written in terms of 
n and a. 

In electromagnetism, having chosen a specific reference frame, one extracts the electric 
field 3-vector from the field tensor F^ by E* = —F t0 , i = 1,2,3, while the magnetic field 

3- vector B' = Fw, where F IW = e Atl/Q ,^F“ /3 /2 is the dual of F^. Given u one can also define 
the two 4-vectors B M = F ,,u u y = q(v • B, B — v x E) and E^ = F^ v u u = y(v • E, E + v x B). 
Note that B* = B l and E* = E l in the reference frame defined by u 11 . In terms of the 

4- vectors B 11 and the field tensor has the decomposition 

= e^ Xa B x u a + (u^E v - u u E») , (9) 

a form valid for any timelike 4-vector uO One can also reverse this process by taking B 
and E M to be fundamental, and then defining the field tensor F^ via (9). In this case, 
different values of B ^ and can give the same field tensor, for one can add any quantity 
proportional to to either 4-vector while leaving the field tensor unchanged; however, if 
the constraints E x u\ = B x u x = 0 are imposed, then the representation is unique. This 
multiplicity of representations of the field tensor will prove important later. 

In MHD one eliminates the electric field from the theory, if necessary using Ohm’s Law 
to express it in terms of the fluid velocity and magnetic field. In a relativistic context, this 
is done by setting = F^ x u x = 0, which gives E + v x B = 0 (i.e. Ohm’s Law) and, in a 
specific reference frame, 

B^ = 7 ■ B, ^ + v (v • B)^ . (10) 

For convenience b^ = B^/^/Att will be used, in which case the MHD field tensor and its dual 
have the forms 

= v / 47t e^ vXlT b\u a and F^ = VA^ (b'\F - u % v ). (11) 

Although (10) satisfies the restriction b x u x = 0, we noted earlier that this condition is 
not needed for a representation of the form of (9). One can, in fact, construct a family of 
vectors 

h? = b^ + au» ( 12 ) 
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where a is an arbitrary scalar field and now, in general, = «^0. The field tensor F 
and its dual F^ u are unchanged when written in terms of hfi, i.e. 


F^jyJ 4^ = e^ x<T b x u a = e^ Xa h x u a (13) 

F^/VAn = \Fu v - u% u = hPu v - u»h u . 


Because only appears in the equations of relativistic MHD via the form (11), one can just 
as easily use the quantity h/fi choosing a in order to give it some useful property. When 
constructing an Eulerian action principle (with covariant Poisson bracket) for relativistic 
MHD it will prove fruitful to do so. The quantity 5^6^, which appears in the stress-energy 
tensor and will be seen in Sec. Ill to appear in the action, evaluates to 


bjf = — (E 2 — B 2 ) 
^ An y ’ 


1 

An 


B • B . 2 

- 2 —I- (v ■ B) 2 

ryA 



2 

rest ? 


where ‘rest’ indicates a rest frame quantity. Thus the 4-vector is spacelike. However, 
since = b^ + a 2 , the status of hfi will depend on a, remaining spacelike for small a. 
Each equation of relativistic MHD can be written as the vanishing of a divergence: 


d^(nu^) = 0 

(14) 

d^au 1 *) = 0 

(15) 

<9„ F ,LU = 0 

(16) 

<9„ = 0. 

(17) 


Equations (14) and (15) express conservation of particles and entropy, respectively. In 
addition, (16) provides the equivalent of the homogeneous Maxwell’s equations; however, 
one cannot call them Maxwell’s equations without qualification, as the constraint F^ v u v = 0 
is already built in when one expresses F^ v in terms of 6^ or h/h 


d v {V x u v - u^b u ) = d u (Wu u - u^h u ) = 0 . 


This expression, of course, is the same whether b M or hfi is used, as the quantity a cancels 
out. Equation (17) gives conservation of stress-energy, where the stress-energy tensor T 7 "' 
is considerably more complex when written in terms of h M rather than b^\ 



where the fluid and field parts are 


7y=(p + p)«V-p 9 '", 

t em = ^ (f“ x F x * + i sTF^fA 

= -Wb v - ( b x b x ) u^u v + ^g^b x b x (19) 

= -h»h u - (h x h x ) u <J u v + [h x u x ) {hV + u»h v ) + (h x h x - (h x u x ) 2 ^j , (20) 

respectively. Equation (19) is obtained by substitution of the first of Eqs. (11) and making 
use of the orthogonality condition b x u x = 0, while (20) follows from (13) without orthogo¬ 
nality. We emphasize that, despite appearances, T^ V M does not depend on one’s choice of a. 
The field part T^ M depends on or hfi only through the tensor F^, in which, as previously 
noted, a cancels out. Lastly, we note it can be shown that this system preserves b^u^ = 0 
and = 1. We next turn to the problem of devising an action principle for this system. 


III. COVARIANT ACTION PRINCIPLE FOR RELATIVISTIC MHD 


The covariant Poisson bracket formalism of Ref. 22 requires two parts: i) an action S that 
is a covariant functional of the field variables and ii) a covariant Poisson bracket { , } defined 
on functionals of the fields. Instead of the usual extremization SS = 0, the theory arises 
from setting { F , S'} = 0 for all functionals F, which is in effect a constrained extremization. 

A general Poisson bracket for fields T has the form 


{F,G} 


" , 5F SG 
Z hT ’ 


where 5F /5T is the functional derivative, dz is an appropriate spacetime measure, and J is 
a cosymplectic operator that provides {F, G} with the properties of antisymmetry and the 
Jacobi identity. Thus 

{F,S} = 0 V F j|| = 0. (21) 

If J is nondegenerate, i.e., has no null space, then (21) is equivalent to JS'/JT = 0 and 
the covariant Poisson bracket formalism reproduces the conventional variational principle. 
However, of interest here are matter models like MHD, which when written in terms of 
Eulerian variables possesses nonstandard or noncanonical Poisson brackets (see e.g. Ref. 23), 
for which J possess degeneracy that is reflected in the existence of so-called Casimirs (see 


9 


Sec. VA). For such systems the covariant Poisson bracket naturally enforces constraints. In 
field theories that describe matter, understanding the null space of J may be a formidable 
exercise 24 , and finding nondegenerate coordinates, which are expected to exist because of 
the Jacobi identity, may only serve to obscure the structure of the theory. 

A variation that preserves the constraints, referred to as a dynamically accessible variation 
in Ref. 25 (see also Ref. 23), can be represented as 

6* da = {*,G}, (22) 

for some functional G, whence 

ss = Jdz^W DA = l'dz^{<t,G} = {S,G} = 0, 

which shows directly how the Poisson bracket effects the constraints without them being 
explicitly known. 

A. Action and functional derivatives 

We construct our action S in a straightforward fashion: 

S[n,cr,u,F] = Jd A x Q(p + p)u x u x + ^(p - p) - -^-F XrT F Xa ^j (23) 

S[n,a,u,b\ = ^ J‘d A x ((p + p - b x b x )u x u x + p - p) (24) 

S[n, a, u, h] — ^ Jd 4 x ((p + p - h a h a )u\u x + (h\u x ) 2 + p - pj . (25) 

Equation (23) is the sum of the fluid action of Ref. 22, where thermodynamic variables p 
and p are considered to be functions of n and a, together with a standard expression for the 
electromagnetic action. 

In (24) the MHD expression of (11) has been substituted into and finally in 

(25) we obtain our desired form in terms of hA Observe that the integrand of (24) when 
evaluated on the constraint U\U X = 1 is the total pressure, fluid plus magnetic, p+ \b\b x \/2. 
This choice of action will be seen to give the desired field equations when inserted into the 
covariant Poisson bracket. 

From the action of (25) one derives a momentum m^ by functional differentiation, 

= — = (p + p - h a h a ) Uf, + (h x u x ) h M + ah M . (26) 
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The quantity 


H = p + p - h x h x 


(27) 


is a modified enthalpy density. If cm M is small compared to b^, /i M will be spacelike, leaving 
fi positive. 

Since u M and b M are independent of a, expressions for them solely in terms of and 
can be obtained. Using a = /iaw a , which follows from (12), and u M = (m M — ah^) //x, which 
follows from (26), we have 


Q? 


= /i a w A = — ( h\m x — ah\h x ) 


Then, solving for a gives 


a = 


h\m > 


(28) 


/i + h a h a 

Equation (28), incidentally, shows that a can be written entirely in terms of the field variables 
and /dk Thus, one can also write the variables 6^ and entirely in terms of the new 
ones: 


= 


m h 


h\nv 


/i + h a h?) 

( h\m x ) 2 




W = /U I 1 + 


h\m y 


(29) 




/i(/i + h a h a ) 2 J n(fi + h a h a ) 

Equations (29) are not invertible. To see this consider a local frame in which v = 0, 
i.e., one where u M = (1,0) and b M = (0, B)/\/47r. In this frame = (a, B/\/47r) and 
= (p + p + 5 2 /4vt, aB/v / 47r). Given any value of cc these equations are compatible with 
(28), but produce the same rest frame values of b M and v?. Thus, Eqs. (29) are not one-one. 
We will explore this degeneracy, which provides a kind of gauge condition, more fully in 
Sec. V. 

Now we are in a position to obtain our action in terms of the variables m^ and h 11 , which, 
due to the form of the upcoming bracket (38), are the appropriate variables for the action 
principle: 

1 (mx* (W) 2 \ (30) 


S{n,o,m,h]=-Jd' x ^ ^ + Khc) 


Upon introducing the “mass” matrix 


M = 


H + aa 


(31) 


a 
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(30) can be written compactly as 


S =-Jd A x (V x - M- 1 - + h x h x - a 2 + p - p) 

— - I d A x ( u x m x + b x h x + h\h x — a 2 + p — p) 


(32) 


= - J d A x (<h A • M -$ A + b x b x +p-p) 

= - f d 4 x (u x m x + b x h x + b x b x + p - p) 


where = ( m x ,h x ), $ A = ( u x ,b x ) and • indicates summation over the 2x2 matrix A4. 
However, because the mass matrix (31) depends on the held variables via /i and a, as given 
by (27) and (28), the expression (30) is superior for calculations; in addition, the mass 
matrix is inconsistent in units, so it would have to be normalized before, say, eigenvalue and 
eigenvector calculations could be done. One possible normalization is given in (57) below. 

After taking variations of the action, one may impose the constraint u x u x = 1. In terms 
of the momentum this constraint becomes 


1 = u x u x 


1 

h 2 


m x m x — 2 


(■ h x m x ) 2 
p + h a h a 


+ 


(h x m x ) 2 
(/i + h a h a ) 2 



(33) 


Thanks to the relations (29) and (33), all functional derivatives of the action of (30) can be 
reduced to simple expressions, provided (33) is applied only after functional differentiation. 
To start with, 


5S _ f m x m x ( h x m x ) 2 ( h x m x y \ dp 1 dp 1 dp 

Sn 2 p 2 + 2p 2 (p + h a h a ) + 2p (p + h a h a ) 2 ) dn + 2dn 2 dn 

dp 
dn 


(34) 


Similarly, 


5S_ _ _dp_ 

da da 


( 35 ) 
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The remaining functional derivatives are 


5S 

5m v 

5S 

W 


m v 


( h x m 


X\ 


H p(p + h T h r ) 

- A (h A m A ) 2 


m x m 


h,. 


hp Uu - 

h„ 




/i 2 (/i + h a h‘ 


(h x m x ) 

/i(/i + /io-h 0 




1 + 2 


(h A m A )' 


(h A m A )' 


+ 2 (/i + h a h a ) (n + h^) 2 

(h A m A ) 2 t ( h x m x ) 


(h T h T ) K 


h,j 


p 2 (p + h a h a ) p(p + h a h a ) 


-m v 


1 + 


(h x m x y 


/i(/i + h a h a ) 2 


hi. 


( h x m > 


p(p + h a h° 




K. 


(36) 


(37) 


The compact result 5S/5h u = b v gives a meaning to h u : it is a conjugate to b u , just as m u 
is to u v . 


B. Covariant Poisson bracket and field equations 


The covariant Poisson bracket for relativistic MHD is obtained by extending the nonrel- 
ativistic bracket of Refs. 20 and 21 to spacetime. This is done by merely summing over the 
four spacetime indices instead of the three spatial ones and altering a few signs. However, 
a difficulty arises in choosing an appropriate equivalent of the nonrelativistic momentum 
and field, because the 4-vectorial equivalents of M = pv and B will no longer produce the 
correct equations. Instead, the 4-vectors m v and h v provide the appropriate replacements, 
giving the relativistic MHD bracket 


{F ' G} yf 


= / drx n 


5F 5G 
5m„ M 5n 


5G a 5F \ 


' 5F 5G 
5m,, M her 


m„ 




5F n 5G 

■On 


5mn M 5m u 


5G „ 5F 

-a, 


5m ^ ^ 5m u 


+ h v 


A 1 A 1 

' 5F 5G 
5m, 11 5h v 5m 


5G 5F 
5m„ M her 


5G n 5F 

~ d “W 


(38) 


d„ 


5F \ 5G 


- d, 


5G \ 5F' 


‘ 5m„ ) 5h u V ^ 5m v ) 5hy 
The bracket is complicated, but one can derive the equations of motion fairly quickly, thanks 
to the simple functional derivatives, as obtained in Eqs. (34), (35), (36), and (37), for the 
action of (30): 


5S__dp 5S__dp 5S_ V 5S 
5n dn ’ her da ’ 5m u ’ 5h v 
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where and b^ here are shorthands for their expressions in terms of the fields and /d 1 
as given by (29). 

Using F — f d 4 x n(x)<5 4 (x — x 0 ) in {F, -S'} = 0 gives, after an integration by parts, 


d^nu^) = 0 , 


which is the continuity equation (14), evaluated implicitly at xo ; however, since that point 
is arbitrary, the result holds for the entire spacetime. Going forward such niceties involving 
delta functions will be skimmed over. In the same manner one also finds the adiabaticity 
equation (15) from a a variation. 

The hh variation gives 


d^hV) - h u d u u» = 0. 


(39) 


The above equations are not Maxwell’s equations, although they are analogous to the non- 
relativistic equation (4), since they correspond to £ u h tl = 0, the Lie-dragging of the four¬ 
dimensional vector density by The theory obtained from the variational principle 
can be viewed as a family of theories, only some of which correspond to physical systems. 
However, if = 0, then one obtains the usual form of relativistic MHD. The situation is 
exactly analogous to that in nonrelativistic Hamiltonian MHD, which can describe systems 
with V'B/0: in both cases, the physical systems are a subset of the full class of systems 
described by the formalism. In the nonrelativistic case the condition V • B = 0 is maintained 
by the dynamics and the similar situation that arises for hh will be shown in Sec. VB. There 
also exists an alternative bracket that builds in = 0, given later in Sec. IV A, where the 
constraint is enforced by the bracket’s Jacobi identity. In any event, with thus specified, 
we can subtract a term u^d u h v from (39), giving the usual equivalent of Maxwell’s equations 


0 = ^(/iV - u^h u ) . 
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Finally, the m x variation gives, after some work, 

0 = — — ad 1 ' + m/ (u u ) + d v (m M u I/ ) 

+ h v & 1 ( b v ) - d v (W) 

= —d^p + (pu u + ( h x u x ) h u ) d^u 1 ' + d v (pu^u 1 ' + ( h\u x ) h^u u ) 

+ hyd 11 (h u - ( h x u x ) u") - d v (h u h IJ - (h x u x ) h v u») 
= du( {p + P~ (h x h x )) uV + ~P + \ (hxh* - (h x u x ) 2 ^j 

- Wh v + ( h x u x ) (Wu v + vPh v ) ) , 




which is the momentum equation (17). Having been derived, it can be replaced with the 
much simpler, equivalent version involving b IJ ’. 

Now we have shown that the covariant Poisson bracket formalism produces the held 
equations of relativistic MHD. In Secs. VB and V C we will probe more deeply the cor¬ 
respondence between the variables (m M , h M ) and (u^,6 M ), exploring in particular how one 
might use the held equations in practice. First, however, we will demonstrate several ways 
in which the bracket formalism can be modified. 

IV. ALTERNATIVE BRACKETS 

In this section we present additional Poisson brackets. The hrst (Sec. IV A) adds an extra 
constraint to (38), the second (Sec. IV B) rewrites the magnetic parts in terms of a tensor 
potential, the third (Sec. IV C) recasts these terms in differential-geometric language, and 
the last (Sec. IV D) incorporates an arbitrary background gravitational held. 

A. Constrained bracket 

Consider the magnetic held part of the bracket of (38), 



(40) 


Just as the nonrelativistic bracket of Ref. 20 and 21 has a counterpart in Ref. 17, the terms 
(40) have an analogous relativistic counterpart that requires divergence-free magnetic helds, 
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i.e. an such that = 0. This relativistic counterpart is simply given by an integration 
by parts of (40) and making use of = 0, i.e., 


{F, G}o h= o 




SF SG 
Srrtfj, ^ Sh v 


5G_ 5F \ 

Sm M Sh v ) 


+ h? 


5F\ SG 
''SlF) Sm v 


' SG\ SF 1 
11 Sh") Sm v 



(41) 


The bracket is identical to (38), but for the swapped functional derivatives in the final line. 
The action (25) is unchanged, as are the n equation (14) and the a equation (15). The h M 
gains an extra term, and may be written directly as the Maxwell-like equation 


d v (h»u v - u»h u ) = d v T ,,v = 0 

without yet imposing a condition on h^. Finally, the equation for ends up with a couple 
fewer terms than before, yielding 

d v V w + (h 11 - {h a u a ) u^)d u h v = 0 , (42) 

where T is the (unchanged) stress-energy tensor (18). 

However, unlike the prior bracket (38), the bracket (41) fails to satisfy the Jacobi identity 
unless the condition d v h v = 0 holds, as is shown in the Appendix. On the plus side, the 
momentum equation (42) is now reduced to its desired conservation form; on the minus side, 
the bracket is defined on a smaller class of functionals than our original bracket (38). The 
original bracket always yields a momentum equation that is not only in conservation form, 
but also independent of a; however, it will yield differing magnetic equations depending on 
a, and only those corresponding to d v h v = 0 produce a Maxwell-like equation. 

We regard the first bracket (38) to be superior, for then relativistic magnetohydrody¬ 
namics may be regarded as a specific example of a broader class of (mostly non-physical) 
dynamical systems, some of which may be of theoretical interest. For instance, in the nonrel- 
ativistic case the broader class have been argued to be superior for computational algorithms 
(see, e.g., Ref. 19), and although similar numerical techniques have been used for numerical 
relativity (e.g., Ref. 5), our formulation provides a fully covariant form analogous to non- 
relativistic MHD that may provide advantages. Moreover, they may correspond to exotic 
theories, such as those including magnetic monopoles. 
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B. Bivector potential 


The divergence-free condition can be made manifest by introducing an antisymmetric 
bivector potential A VIX such that 

h " = d v A v ». ( 43 ) 

Such a representation gives rise to a gauge condition — > A ,w + for chosen -0^; 

such gauging could be useful, but we will not explore this further here. 

Assuming F[h\ = F[A], i.e. functionals of the bivector potential obtain their dependence 
through h , we obtain 

5F = d A x—5h» = / d 4 x — — 8A Ufl = 5F. ( 44 ) 

/ 5h^ J 5A"v v ' 


Relate 5h M to via (43) and insert h/A = d v 5A u ^ into the second equation of (44). Even 

assuming 5A UM is arbitrary, it only picks out the antisymmetric part of what it is contracted 
with, so we obtain the functional chain rule relation 


5F 1 f SF 8F\ 

5A^ 2 \ u 5h» ^5h")' 

Inserting (45) into (41) gives the compact expression 


{F,G} a -.= 2 d 4 x (d a A°n 


5F 5G _ 5G_ 8F \ 
5m M 5A u i x 5m M 8A V ^ ) 


We will use this form in Sec. V A, where we discuss Casimir invariants. 


(45) 


(46) 


C. 3-Form bracket 


For nonrclativistic MHD we observed in Sec. IIA that the magnetic equation may be 
written <9B/<9 1 + £ V B = 0, where £ V B is the Lie derivative of the vector density B dual 
to a 2-form. Thus one can write B l = e l Fu}j k and c Oj k = B l 6ij k / 2, where i,j, k = 1,2,3. 
In terms of the 2-form the equation becomes dco/dt + £ v oj = 0, with £ v now being the 
appropriate expression for the Lie derivative of a 2-form in three dimensions (e.g., Ref. 26). 
In n dimensions, an (n— l)-form has n independent components. This suggests we can 
introduce the dual 3-form for relativistic MHD as follows: 


(x(y h and 


h s = -e a ™u afh . 
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which shows that h M is a vector density because it is the contraction of the tensorial three- 
form with e a/ 3 7 5 a relative tensor of unit weight. From the above it follows that the 3-form 
equation of motion is given by du/dt + £ u ui = 0. If we denote by F£ the 4-vector given 
by 5F/5m then the magnetic portion of the Poisson bracket in terms of the 3-form can be 
compactly written as follows: 

Although similar expressions in terms of Lie derivatives exist for all terms of all brackets, 
we are concentrating on the magnetic terms, which written out are 

Gm G^d/jOJup 7 + UJ^p^daG^ + OJa^dpG 1 ^ + OJap^d^G^ . 

The transformation from the bracket {F, G}h of (40) to that of (47) follows from a chain 
rule calculation similar to that described in Sec. IV B. Thus, it satisfies the Jacobi identity 
because {F, G}h does, as shown directly in Appendix A. 

Relativistic MHD has a natural 3-form dual to viz. F\ a u v + F au u\ + F u \u a , which 
follows from the definition b 1 ' = - v /47re / “ /Acr F\ a u v /2 with = 0 and F^u" = 0. The 3-form 
dual to h M can similarly be represented as <jJ\ av = [F\ a w v + F au w\ + F v \w a ) /6, where 
W/j, = {lFu^-ah^) / (b\b x ) is designed so that = 0 and w^ = 1 and evidently u>\ au h^ = 

0. Observe can be written in various ways using (29), (28), and other expressions. 

The Jacobi identity for the bracket with (47) does not require closure of the 3-form. 
However, if the 3-form u is exact then it can be written as the exterior derivative of a 2-form 
as follows: 

9 a Af 3 7 T dpA'ya T bjyAop 

and one can rewrite the bracket in terms of A /ll/ . Instead of writing this out, we observe the 
bivector potential is given by 

A v>1 = 1 e^Aar 

and so the closed 3-form bracket is essentially given by (46). 

When the 3-form co a p 1 is exact we have, for any 3-surface H in our four-dimensional 
Minkowski space-time, Stokes’ theorem 



18 




where f Q u contains the notion of ‘flux’ in this setting. If 12 contains a time-like direction, 
we can write this as a conservation law, but such 3 + 1 splittings will not be considered here; 
instead, we refer to Ref. 22. 


D. Background gravity 

Now we generalize the full formalism to curved spacetimes. In this context, the equations 
(14) - (17) are now written 


(nu^).^ = 0 

(48) 

(au^).^ = 0 

(49) 

F' w . u = 0 

(50) 

T^. u = 0, 

(51) 


where the denotes covariant derivative. 

Three modifications to the previous action principle are required: (1) because all inte¬ 
grations have tensorial integrands, the integrations must take place over a proper volume 
y/—g d A x ; (2) h M should be treated as a contravariant vector, and as a covariant one, 
befitting their definitions (note that treating them any other way would introduce extra 
factors of g^ v into the bracket); (3) functional derivatives should be defined in a way that 
makes them tensorial. Specifically, for a held variable v, one implicitly defines the functional 
derivative via 


d_ 

de 


e=0 


f 5F 

F(v + e8v) = d A x —8v yf—g . 
J OV 


The action is now 


1 /' , 4 ( g Xa m\m a 

b — — I d x 

F 


(h x m x ) 2 

/i(/i + g\ a h x h a ) p P J 9 


and its functional derivatives are 


dS___df^. _5S__ fj, . $S__ 

Sn dn ’ 5a da 1 ’ 5hv ' 9iV 
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Finally, the bracket becomes 


{F,G} = fd A xV=g( 


( ( 5F „ 5G 5G 5F\ f 5F „ SG 5G „ SF\ 
U 'Sm^ln ) +a ) 


m„ 




SF 0 SG 

■a a 


SG 0 hF 

o, 


5mp ^ 5m u 5m ^ M 5m, 


, ^ a ^ ^ a ^ 

1 5m° il 5h v 5m 5h v 


( 52 ) 


a 


hF \ hG 


/ 5h u 


5, 


SG \ 5F 


^ 5m,, ) Sh v 


The d M operators inside the bracket are still just partial derivatives, but the presence of the 
metric will tend to convert them into covariant derivatives (see e.g. Ref. 27, Ch. 21). After 
an integration by parts, the variation {F, S'} = 0 of the test function F — f d 4 xn(x ) <5 4 (x — 
%o) S ives 

(nu^^g) = v 73 g(d M ( nu M ) + nu u T^J = (n/) ;;j = 0 , 

with a similar result obtaining for the a variation. The h M variation once again requires 
special attention, as it gives 


d u (h?u v y/=g) - h v (d v ,F) ^Tg=^r g (hV u + h\u v ~ h v u>' v + h»u x T\ u ) = 0 . 


This time we choose a so that hJ 1 .^ = dpJrF + h^F^ = 0. Similar considerations apply 
to this choice as in the special relativistic case. Subtracting this expression and combining 
like terms then gives, with = h^u v — h u u 


d u F' JV + F^F\ V + F uX F^ uX = 


F ,,lJ 




0. 


Note that the third term is zero by the antisymmetry of T^ and the symmetry of the 
covariant indices of F' J uX . 

Finally, one obtains the momentum equation (51) by varying the test function F = 
Jd A xg^ u m v 5 A (x — Xo) x /—g. This derivation is lengthy, and will only be summarized here: 
(1) the partial derivative terms appear, and combine, exactly as in the special-relativistic 
case; (2) the T ,jX F u Xv terms come from taking the partial derivatives of \/—g\ (3) the T uX F^ uX 
terms come from derivatives of extra factors of the metric g^, some of which come from its 
inclusion in the test function, others of which come from 5S/5h M = g^Jf. 

We conclude with an important note. While we constructed the above formalism to han¬ 
dle curved spacetimes, it also applies to flat spacetimes with arbitrary coordinate systems, 
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such as cylindrical, spherical, or toroidal coordinates. The nonrelativistic version may be 
generalized the same way (altering volumes d 3 x to proper volumes y/gd 3 x), thus solving the 
problem of MHD coordinate changes in a pleasantly general way. 

V. DEGENERACY AND SETTING THE GAUGE 

Now we consider various issues pertaining to degeneracy. In Sec. V A we obtain Casimir 
invariants, showing that the action S is not unique. Then in Sec. VB we further explore 
the noninvertiblily of the transformations from (u M , 6 M ) to (m M , /ffi). Finally, in Sec. V C we 
discuss the how the divergence-free condition on h M can be constructed for any problem. 

A. Casimirs and degeneracy 

As noted in Sec. Ill the covariant Poisson bracket possesses degeneracy and associated 
Casimirs. A functional C is a Casimir if it satisfies 

{F,C} = 0 V F. (53) 

Equation (53) should not be confused with the variational principle of (21), {F, S'} = 0 for 
all functionals F. for the former is an aspect of the bracket alone, and provides no equations 
of motion. Because of the definition of C, the action S is not unique and can be replaced 
by S' + A C for any Casimir C and any dimensionally appropriate number A. 

Turning to the task of finding Casimirs, we use (53) to provide functional equations for 
the Casimirs. Although difficult to solve in general, some explicit solutions can be found, 
facilitated by our knowledge of Casimirs for nonrelativistic MHD 21,28 . First, it is easy to 
obtain a family of what we call the entropy Casimirs, 

Cs = Jd i xnf(a/n) , (54) 

where / is an arbitrary function. In the nonrelativistic case this is a generalization of the total 
entropy, for if / = <r/n and a is the entropy per unit volume then f d 3 x nf(a/n) = f d 3 x a 
is the total integrated entropy. 

Next we seek a Casimir that is a relativistic version of the cross hclicity J d 3 x \ ■ B. For 
nonrelativistic MHD invariance of cross helicity requires a barotropic equation of state and 
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V-B = 0, so we make analogous assumptions here. We assume p has no dependence on a, 
and we implement the analogue of V • B = 0 by using the antisymmetric bivector potential 
of Sec. IVB, h** = <9 7 A 7/ b ensuring that d^hJ 1 ' = 0. Using the bracket of (46) it is easy to 
show that the following generalization of the cross helicity is a Casimir: 



(55) 


This quantity ceases to be a Casimir when the divergence h M is nonzero. Observe that on 
the constraint U\U X = 1, the integrand of (55) can be written as /n = m fl h ll /n = 

a(p + p)/n, which follows from (28). Since a does not exist in the original (w M ,6 M ) theory, 
this Casimir is a quantity tied to the covariant bracket theory in terms of (m M , hk). 

One also expects the existence of a magnetic helicity Casimir, but the nature of linking 
in four dimensions makes the situation complicated. Relativistic generalizations of magnetic 
helicity have been found in Refs. 16 and 29, but we have yet to demonstrate that a quantity 
like either of these is in fact a Casimir. We also anticipate the existence of additional 
Casimirs that are generalizations of the nonrelativistic ones found in Refs. 30 and 31, but a 
full discussion of Casimirs will await a future publication. In any event, because of the form 
^da as given by (22), we can be assured that the extremization of our covariant bracket 
variation preserves any Casimirs that exist. 

B. Gauge degeneracy 

In Sec. Ill A we noted that Eqs. (29) are not invertible. This lack of invertibility, which 
arises from the gauge freedom associated with a, can be understood in greater generality. 

Because the degeneracy is not associated with the thermodynamic variables p and p, we 
move them out by introducing the following scaled variables: 

h = (VP + p) h , m = (p + p)m , b = (y 'p + p)b, u = u, a = (^/p + p) a , 

In terms of these variables (29) becomes 


$ = wr 1 . t, 


(56) 


with 



b 1 —a p + a 2 


—a 


u V ot 2 a 

M = 

a 1 


(57) 
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and $ = (u, b) , 4/ = (m, h). The quantity /x = 1 — h 2 is a normalized /x, and the quantity 
d satishes a = fh v h v = u u h u . Varying (56) gives 

- , - dM - 1 - dM- 1 - 

6$ = M- 1 -6V + - • T hd + ——— • T h/x. 

(7Q a/x 


Degeneracy follows if we can hnd a nonzero hT giving hd> = 0. Such would be given by 


dM- 1 - - dM- 1 - 

SV=~M- -VSa-M- • T h/x 


9a 


9/x 


- dM- 1 - - dM- 1 - - 

= -M - - M -<f>hd- M • -4- M-^Sfx 


da 

dM = . =. 

<f> ha H—• $ h/x 


9/x 


da 
= ha 


9/x 


2a 

1 


1 

0 

1 

0 

• $ + h/x 

0 

0 


• $. 


(58) 


Thus from (58), hm^ = ( 2au u + b u )5a + u u 5fi and 8h u = u u ha. Using h/x = —2 h u 5h u = 
—2h u u u Sa = —2aha, the two conditions imposed by (58) are 


5h v = u v hd and hm" = b v hd , 


(59) 


reiterating our earlier point that a can vary while leaving and b M unchanged. 

In terms of the scaled variables the action becomes 

S[n , a, fh, h] = ^ J d A x (m\fh x — (h A m A ) 2 j + p — pj • (60) 

Now if we consider variation of the integrand of (60) with variations given by (59), and 
restrict to the constraint = 1 as given by the scaled version of (33), then the action is 
easily seen to be invariant. Using the scaled action in the form of (32), the integrand becomes 
upon variation (p + p)(u\5fn x + b\5h x ) + h\5h x — dhd, which vanishes upon insertion of (59), 
with the first two terms vanishing individually because u\b x = 0. Thus, degeneracy appears 
as one transitions from (24) to (25). We add that in scaled variables J- ~ u fJ, b l/ — b^u 1 ' ~ 
u tl h l/ — hf l u v \ thus, at fixed u8T ~ u^5h u — 5h tl u l/ = 0. 

However, we also require that the equations of motion in terms of (xx M , V 1 ) stay unaffected 
by the degeneracy in (m M , h M ). This requires d^M = 0, as we earlier discussed in the context 
of the magnetic equation (39). Written in full, this condition becomes 

df, (axx M ) = -d^U 1 . (61) 
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As usual in the case of extra degeneracy, the system now possesses an additional symmetry, 
for one can add to a any solution A a of the continuity equation (A = 0 while leaving 

the dynamics unchanged. This is not as powerful as choosing a freely, as (59) seemed to 
imply, but we will show in Sec. V C that it is nearly as powerful. We hope to further explore 
the consequences of this new symmetry in future work. 

Our system’s degeneracy is related to the adaptation of Goldstone’s theorem 32-35 de¬ 
scribed in Ref. 36, where it was proven in the context of degenerate Poisson brackets with 
Casimir invariants that nonrelativistic Alfven waves associated with degeneracy can be 
thought of as an analog of Goldstone modes. A similar interpretation arises here in this 
covariant relativistic MHD setting, but discussion is beyond the scope of the present work. 

C. Setting the gauge 

Given a relativistic MHD problem posed in terms of (iG,^), we must determine the 
associated problem in terms of and this requires the determination of a, which 

amounts to setting the gauge so that d^h^ = 0. Doing so may seem difficult on first sight, 
but in fact turns out to be simple. Since this idea sits at the crux of our formalism, we will 
explain it is some detail. 

Posing a relativistic MHD problem requires one specify (u M , 6^) as well as n and a on 
a space-like 3-volume, Q G T>, where is T> is our four-dimensional space-time. In addition, 
a physical problem will have initial conditions that satisfy U\U X = 1 and U\b x = 0. Using 
u°d a = djdr where r is the proper time measured by an observer comoving with a flow line, 
one can choose r = 0 to correspond to the state specified on and then propagate values 
off of by using the equations of motion to determine db^/dr, du^/dr, dn/dr , and da/dr 
at r = 0. This is the standard scenario for a Cauchy problem, and many references for both 
MHD and relativistic fluids (e.g., Refs. 2 and 4) describe this in detail. One can imagine 
an exotic flow in which there exist spacetime points not connected to D by any flow lines; 
however, a modest boundedness condition excludes such cases. 

The present situation is complicated by the fact that given on Q at t — 0 we must 
also have that = 0 for all time, in order for our (m /1 , h dynamics to coincide with 
the physical (tC, b^) dynamics. Fortunately, d^h^ = 0 is maintained in time if it is initially 
true on D. To see this we act on (39) with d M and obtain d u {u u d^h^) = u u d^d^) + 


24 


(d u u u )(d lJ ,h ,Jj ) = 0 or 


^P- + (9u u‘’)(d„h‘‘)=0, (62) 

an equation analogous to (6) for nonrelativistic MHD. From (62), one concludes that if 
d^Jn/ 1 = 0 on O at r = 0, then d M hd remains zero for all time. Thus, one can solve the 
(m M , h^) equations and uniquely obtain the (w^, b via (29) - provided one can ‘set the 
gauge’, i.e., find an a such that d M hd = 0 on hi at r = 0 consistent with the (id 1 , n, a) of 
our posed problem. 

We will first consider a special example of setting the gauge, corresponding to the case 
described in Sec. Ill A. We are given the MHD problem with initial conditions v(0,x) = 0, 
i.e., w^Ojx) = (1,0) and fd(0,x) = (0, B(0, x))/\/47r on the space-like 3-volume D with 
coordinates x, and we wish to obtain an hd(0, x) = (a, B/V^7r) and m M (0,x) = (p + p + 
B 2 /Att, aB/\/47r) such that <9 M hd( 0, x) = 0. Denoting do a = a t , etc., gives the condition 

0 = <9 M hd(0, x) = —= ( 7 tv ■ B + yvj • B + qv ■ B t + cq-y/Urj + V • h , (63) 

where h is the spatial part of h A Evaluating (63) on the initial condition gives 

0 = Vf • B(0, x) + o/.t (0, x)\/47r + V • B(0, x), 


whence, with V • B(0,x) = 0, we conclude that 


0 = v t • B(0,x) + cq(0,x)\/47r = --'S/p - B(0, x) + a t ( 0,x)\/47r, 

P 

using the MHD momentum equation in the last step. Thus a t ( 0, x) = (a/Ht p)~ 1 S7p ■ B(0, x) 
on fl will assure c^/d 4 = 0 for all time. Observe, a(0, x) has not been specihed - we are free 
to choose it as we please. In doing so we will obtain different initial conditions m tt ( 0, x) and 
h /4 (0,x) and these can be chosen for convenience. Finally, if we solve our equations for md 
and hd and obtain their values at any later time, insert them into (29), then values of 
and fd thus obtained are solutions of the relativistic MHD equations. 

Now let us consider the general case, beginning with the expression 


= djf + ^(rnd 4 ) = d f Jf + d^anu^/n) = d^ + n 0) 

where the last equality follows from (14). Upon contracting d p (b fl u l/ — u^b u ) = 
obtain 


d v b v = u v 


db v 


-b v 


du v 

lh' 


(64) 

0 with un we 

(65) 
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Consequently, (64) and (65) imply 


d /a 
dr \n 


b v du v 
n dr ’ 


a{r f ) = n(r f ) 


rTf b u du, 


n dr 


■dr + a(0), 


( 66 ) 


where the above requires one integration per flow line. 

Thus the freedom in a reduces to a choice of a on the initial surface Q, its value at any 
later time being found by solving the Cauchy problem. Furthermore, even this initial step 
may be rendered trivial. While discussing the condition (61), we pointed out that one can 
add to a any quantity Aa obeying <9 M (Aa w M ) = 0. Reiterating the argument that led to 
(66), we find this becomes 

V 


or \ n ) 


which says that we may choose Aa freely on 0, and the ratio of it over the number density 
will remain constant along flow lines. Given this freedom, why not simply pick Aa = —a on 
the initial surface? So the new a is zero on D, the initial conditions are simply (m M , hP') = 
({p + P + |6 2 |)m m , b^), and a develops along flow lines according to (66). Said integral never 
actually has to be evaluated, for if one solves the Cauchy problem for and (whose 
equations of motion incorporate the condition d^h^ = 0), one can then calculate a via 
(28). Nonetheless, (66) may be useful as a consistency check on calculations or simulations. 
Similarly, the two constraints u^b^ = 0 and = 1 propagate along the flow lines and 
do not need to be enforced explicitly provided they are true on initially, though they too 
remain useful as consistency checks. 

We close this discussion by considering a point that may cause confusion. Given (mF, h M ) 
on we can certainly calculate V • h, and dh°/dr will be determined by the equations of 
motion for (m M , h fi ). Thus, one may wonder how we are free to chose a and da/dr to make 
dpJP 1 = 0. The answer lies in the fact that the (mfl, h M ) system has a solution space that 
includes solutions that are not relativistic MHD solutions, and our procedure for picking the 
quantity selects out those that do indeed correspond - for these the two ways of determining 
d^hP are equivalent. 


VI. SUMMARY 

We have successfully cast relativistic MHD into a covariant action formalism using a 
noncanonical bracket. Along the way, we had to develop a few new ideas with possible con- 
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sequences beyond our current domain: a modified enthalpy density containing a magnetic 
“pressure”, a canonical momentum differing from the kinetic momentum by a magnetic 
term, and a divergenceless magnetic 4-vector possessing a new degeneracy and symmetry. 
We presented several closely related additional brackets, and carefully investigated the non- 
invertibility of the transformations between our original Eulerian quantities and their con¬ 
jugate momenta. Many consequences of our formalism were investigated, but many more 
remain to be covered: for instance, 3+1 reductions, additional Casimirs, the relation to 
Lagrangian action principles, brackets in systems possessing extra symmetry (e.g. spheri¬ 
cal or toroidal), applications to the Godunov numerical scheme in relativity, and conserved 
quantities related to the a symmetry, ft may be objected that we have, as yet, produced 
no practical application for our formalism, though it certainly does possess a certain beauty 
of its own. However, while practicality usually precedes beauty in physics, the opposite is 
sometimes the case, reason enough not to disregard that beauty. 
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Appendix A: Direct proof of the Jacobi identity 

The brackets of (38) and (41) are direct generalizations of the Lie-Poisson form given in 
Refs. 17, 20, and 21 for nonrelativistic MHD, so the Jacobi identity follows from general Lie 
algebraic and functional derivative properties (see e.g., Refs. 21, 23, 37, and 38). However, 
since these will not be known to most readers we include a direct proof in this appendix. 
The Jacobi identity is 

{{F, G}, H} + {{G, #}, F} + {{H, F}, G} = 0 (Al) 
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for the two brackets (38) and (41). 

When expanding the expression (Al), many terms will contain second functional deriva¬ 
tives, for instance 


nh } 


5G 
5m,, 


d„ 


5 2 F 


5h x 5m^J 5n 


d, 


5H 


Thankfully, by a theorem in Ref. 21, all such terms cancel for any antisymmetric bracket. 
Thus we only have to worry about those terms containing only first functional derivatives. 
Starting with the bracket (38), the needed terms are thus 


HF,G} 

5n 

H A, G} 
5a 

5{F,G} 
5mn 
5{F,G} 
5h^ 


5F 5G 
5mn M 5n 
5F 5G 
5m „ M 5a 
5F 5G 

Ou 


5G 5F_ 
5m M M 5n 
5G_ 5£ 

5m,, ^ 5a + 


5m, 

5F 


d, 


5m^ 
5G 


5G 

5m u 

5G 


, SF 

O v - -h . . 

5m,, 


(A2) 


* SF 
d u —— + 


5m u 5hv 5m u 5h f 

„ 5F 5G „ 5G 5F 

d„ -—— - d „-—— + 


' ,J ' 5m u 5h v 11 5m u 5h u 
with similar expressions for the other two permutations of F, G, and H. Beginning with this 
expression, it is to be understood that, in the absence of parentheses, the gradient operators 
act only on the term immediately to their right; when they are followed by an expression 
in parentheses, they act as normal. This convention will remove many superfluous symbols. 
The ellipses at the end of each line indicate the terms that may be disregarded thanks to 
the aforementioned theorem. Upon inserting the expressions (A2) into the Jacobi identity 
(Al), all pertinent terms will be linear in the held variables. Each of these four sets of terms 
(one for each held variable) must vanish separately. 

The terms linear in n are: 



( 5F 5G_ 
\5m u v 5m^ 


5G 5F\ 5H_5H (5F 5G_5G 5F\ 

5m u v 5mf,) ^ 5n 5m^ M \ 5m u v 5n 5m v u 5n ) 


+ O 

F,G,H 


(A3) 

where the circle symbol indicates permutation in F, G, and H. Inside the square braces, 

the collected second derivative terms are 

5H 5F 2 5G_ FH_$G_ d 2 $£ _ _^LJ^_g2 5 Ji 

5mn 5m u lll/ 5n + 5m M 5m u ,M/ 5n 5m M 5m u ^ 5n 

5F 5H 2 5G _ $G_FH_ d 2 SF_ 2 

5mn 5m u ^ 5n 5m^ 5m u 5n 5m^ 5m u ^ 5n 
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which vanish due to the fact that second (partial) derivatives commute. The remaining 
terms linear in n, keeping the same order they have in the Jacobi identity, follow: 


5F „ 5G „ 5H 


8m,, v 5m,, '' 5n 


-d, 


-d, 


-d„ — 


8m,, v 5m„ ^ 8n 


5G „ 5F „ 5H 


8m,, u 5m,, M 5n 


-d, 


8H „ JF „ JG 


-d, 


d u -z — + 


5G „ JF „ JF^ JF „ JG „ JF 


® 


d„ 


-d„ — 


5m,, ^ 8m,, u 8n 


5F „ JG „ 5H 


8m,, v 5m,, ^ 8n 


-d, 


8m,, M 8m,, v 5n 


5H „ JG „ JF 


© 


-9, 


9 d~z — + 


® - + 


5m„ v 5n 


8F „ JF „ JG 


-9, 


d v —— + 


5m,, ^ 8m,, u 5n 


5H 5F 5G_ [ 
8m,, v 8m„ M 5n 


5F 5H 5G ( 
5m,, u 5m„ M 8n 


5G_ JF 5_F_ k 
5m„ ^ 5m,, u 5n 


5G 5F 0 5JF 
5m,, 11 5m,, v 5n 


They vanish in pairs, as labeled by the circled numbers. 

So all the terms linear in n have vanished from the Jacobi identity. However, the terms 
linear in a are identical, but with functional derivatives 5/5n replaced by 8/8a. So the a 
terms vanish by an identical calculation. Moreover, the m\ terms do as well: the 5/5n are 
replaced with 8/5m\, contracted with the remaining m\ term outside the square brackets 
of its version of (A3), and the calculation proceeds as before. 

The only terms remaining to be checked are those linear in h x ] unfortunately, there are 
quite a few: 

5H ® 


d x h 


5F 5G 
5m^ M 5m u 


5G 5F 
5mp M 8m u 


d, 


5H 

8m u 

+d x 


d,. 


5F 0 5G 

■a„ 


5G 


5m^^5h x 

5F n 5G 
d„ 


d„ 


5h x 

dF \® 


^ dh x J 


5H 

5m,, 


d„ d. 


5F 5G 
' 5m^ 5hv 


dy 


8mn ^ 5m 


5G „ 5F 

a, 


5m ^ M 8m, 


5H 

w 


-d> 


5H 
5m,, 


5F 5G 
5m„ 5h v 


5G 5F „ 5F 5G „ 5G 5F 

- d - +d - d - 

5m,, M 5h u u 5m„ 5h^ u 5m„ 5h^ 


5G 5F 
5m ^ 5h ft 


+ O 

F,G,H 


The terms labelled by a circled “one” produce a calculation identical to that already per¬ 
formed, and thus cancel. From the remaining terms, we first gather all the second derivative 
ones inside the square braces: 


5H 5G 
5m u 5hv 

5F 5H 
5m u 5h^ 

5G 5F 


d' 


\v 


d’ 


\v 


dl 


5F 
5m ^ 

5G 

5m f, 

5H 




8m,, 5h^ Xu 5m 


5H 5F 5G 

ai 


5m,, 5 Xu 5m 


5F 5G 


dl 


8H 


5m v 5h^ Xu 5m^ 

5G 5H 5F 
oi 


5m,, 8 /F Xv 5m 


5F JF 2 5G 
a 


® 


5m fl 5h l/ Xtl 5m 


5G 5F 


dl 


5H 


8m^5h v Xtl 8m h 


5H 5G 


dl 


5F 


5m ll 5h u x ^8m v 


5G 5H 2 5F 
di 


5m^5h u A/i 5m h 

5H 5F 5G 
a 




5m^5h u A/i 5m 


5F 5G 


dl 


5H 


5m ll 5h v Xfl 8m h 
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They cancel in pairs. Finally, the remaining terms, in the same order and bearing the same 
indices as in the Jacobi identity, are: 


5H 0 5F 0 5G ( 
dx 5^ d "5h^ 


5m h 

5F 


SF[ 5G n 5F® 
+ — — d x - —<9,, 


a SH_ 1 5G _® 

5m„ A 5m,, M 5h v 5m 


5m, 

5G 


' 5m„ "5h^ 

x 6H x 6F{ 
0\- - On 


5H 0 5F 0 5G® 
o\- — d„ 


5by 

5G 


‘ 5mn M 5m v 


5_H_ 5G_ 5F_ 

5h v x 5m„ M 5m,, 


5F 5G 5Ff [ 

dx J^ du Ih^ 


5m 


5F 

+ 7 — d\ 
dm,. 


5m v 

5H 


5h h 


SH 6F® SF 

0\-z -°V7- + TTT^A 


'V 

5H 


^ 5G® 
5rn~ dv 5hi I 


5hv 5m u Srrif, 

5F „ 5G . hi/ 
+ —<9 A ^— 


<9 — 

J/d 1 v 5m„ 


o 


5G 5F 5H ® <WJ 5F 5G' 
5m^ X 5m u 11 5h v + 5m^ X 5m u ^ 5h v 

5G . JiJ . 5F F 

0\~ - 0,y 


5h u 

5H 


5mn ^ 5m v 
~ 5F 5G ( 

o\- — o,j 


5m 

5H 


' 5m„ u 5h^ 


5G 5F 5H 
+ ~ — d x - — d, 


_ a 5G_ d 5F_ 

5m,, X 5m,, M 5h v 




+ 


5m, 

5F 


Sm^ 5hv 

x SG x SR{ 
dx — — d„ 


5h^ 5m„ 5m M 

L 5G. 5H 5F © 
+ -dx- — d„ 


5F 5H 5G_ 
5h u X 5m /1 M 5m v 

,5G 5F 5H 

+ Ti-dx- —<9„ 


5h l 

5F 


' 5mn M 5m v 


5hv 5m v 5m M 

SG 5F^ d 5H_ 
Sh" X 5m„ M 5m,, 


o 


5m,, 5m,, 5h v 


„ 5G 0 5E^> 5H 5G ^ 5F & 
dx t— d v —— + TT-dx- — d u 


5hv 5m,, 5m, 


5h^ 5m,, 5m, 


They also cancel in pairs, establishing the Jacobi identity. This derivation is also valid in 
curved spacetimes, for the functional derivative cancels out a factor of yf—g, and there is no 
integration by parts to catch another such factor. 

Next we will perform a similar calculation for the alternative bracket (41). While the 
same kinds of terms appear as above, there is no longer a complete cancellation. Most of 
the functional derivatives (A2) are unchanged, the only differing one being 


5{F, G} 


5F n DG 

-d,j 


5G _ dF 0 5F 5G 
— du^~ + <9 M - 


n 5G 5F 
d ,—-— + ... 


5h^ 5m v dh^ 5m v dh^ ^ 5h u 5m u ^ 5h v 5m v 

with the ellipsis again indicating terms with second functional derivatives, all of which can 
be disregarded. 

The terms of the Jacobi identity once more appear in four sets, each linear in one of the 
field variables. The n, a, and m x terms involve no derivatives with respect to h x , and are 
thus unchanged: they cancel as before. Only the h x terms differ. They read: 


5H 

5m,. 


d v 


5F 0 5G 
d„ 


5m^^5h x 


+dx 


5G o 
5m 11 
5G 


5F _ 
5m„ ^ 5h v 


d 4 x h x 

5F\ 

sJf) 

5G 


5F 0 5G 

■du 


5m v 5m„ 5m v 5m, 


5G 5F \ n 5H® 
-du- — o„ 


© 


5H 

5m,, 


du dx 


5F 5G 


5hv 5m,, 


d: 


5h x 


5G 5F 


5hv 5m 


d,,_ + du^-— - d, 


5G 5F 


5mn 11 5h u 


5hv 5m, 


5h^ 5m^J 5m v 


5H 


-< 9 ; 


5H 


5F 5G_ 
5m„ ^ 5m,, 


5G 5F 
5m,, M 5m,, 


+ o 

F,G,H 
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The terms labelled with a circled “one” cancel as in the previous bracket. The collected 
second derivative terms are 

_ SH SG 2 SF ® 5H SF 2 SG ® SF 5H 2 6G 

Sm u Srri/j, vX Sh^ Sm u Sm M uX Sh^ hm M Sm u A/i Sh v 

5G 5H 2 hF hi/ hG 2 hF ® hi/ hF 2 hG ® 

Sm^ 8m u X>1 8h u + 8m u Sm^ Xu Sh^ Sm u Sm^ Xu Sh^ + f,g,h 

6F SH 2 SG SG SH 2 hF SG SF 2 hi/ 

“h^h^G A/i h/F ~~ h^h^G V h/F + hm„ 

5H 6F 2 hG FF ^ 2 

hm M Sm v A/i Sh u hm M hm^ A/x hiz^ hm M hm^ AM h/F 

Six terms do not cancel. The other terms (i.e. those that are not second derivatives) are 


SH SF SG ® SH SG ^ SF 
o\—-d u -— + -— d\-^d v 


Sm 

SH 

Sm,, 


Sh 11 hm M 

o 5G o SF 
dx 8m® u 8h^ 


Sm v Sh^ Sm h 

SH _ SF SG 
~ d ^ d ^ ~ 


_ SH_ Sh? SG_ 
Sm v A Sm,, M Sh 1 ' 


SH SG SF_ 
Sm u X Sm lx 11 Sh v 


SF SH SG ® SG SH SF 
Sm,, X Sh u M Sm,, Sm,, X 5h u M Sm,. 


SF SG „ hii 

d\-Fr-o„ 


Sm, 

SF 


Sm u 

SG 


d. 


' Shv Sm 
SH 


-d,j 


SG 


Sm u Sm M 

SF SH SG ® hF hG hF 
hm„ A h/F v Sm„ + Sm,, X Sm„ M h/F 


Smn Shv 


SF SG „ SH 

-d \-— 


SH 


A 4 

hF 

hm,. 


o 5F hG 
dx Sm^Sh" 


d\^rr-d, y 


Sm u Sh^ Sm M 

t SG ^ SF ^SH 
+- —<9 a v—a, 


hm„ hm^ h/F 

SG n SF n SH 
~d\— — u„ 


SG Sh? _ 

Sm„ X Sh v ^ Sm 


SH SF SG 
Sm„ X Sh v M 5m v 


Sm h 


Sm M h/F 


hm„ h/F hm M 
hG hF hF 

-OA“— O v 


SG SH SF 
hm„ X Sm„ 11 Sh u 


Sm,, Sm,, Sh^ 


SH hG SF_ 
Sm„ X Sh u 11 Sm,, 


F 

SG 

Sm, 


d HL d 5 2L 

° X Sm^ Sh v 


SF SG_ SH_ 
Sm„ X Sh v M Sm,, 


°/J, W 1 0 VIILfj, KJIV KJULjy 

This time twelve terms do not cancel. All told, eighteen terms remain, which collect in 
groups of three. Each group reduces to a gradient with a B\ pulled outside the expression. 
The whole Jacobi identity simplifies to 


{{F, G}, F} + {{G, F}, F} + {{F, F}, G} 

SF^SG_ Sh^_SG_SF_ SH_ SG SH SF 
Sm v hm M M Sh v Sm u Sm M 11 Sh v hm^ hm^ M h/F 

SH SG SF SH SF SG _ SF SH SG\ 

Sm u Smn /J Sh u hm„ hm M ^ Sh v Sm u Sm M M 5h u ) 

An integration by parts shows that the Jacobi identity is satisfied if h u u = 0. In a 

curved spacetime, the above expression is the same, except that <Z 4 x becomes ^J—gd A x. 
The integration by parts catches this extra factor, yielding (h v ^J—g\ v = h v , v = 0 as a 
requirement for the Jacobi identity. 
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